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What can be inferred from surrogate data test- 
ing? 

In a recent Letter Palus and Novotna ||l[ reported sta- 
tistical evidence based on surrogate data testing for lin- 
earity that a driven nonlinear oscillator is the mechanism 
underlying the sunspot cycle. 

While their result might be true we doubt the formal 
correctness of their conclusion. Surrogate data testing for 
linearity [^ tests the null hypothesis that a linear, Gaus- 
sian, stationary, stochastic dynamical process underlies 
the data, including a possible invertible, static nonlinear 
observation function. To perform the test a feature is 
chosen that can capture a violation of the null hypothe- 
sis. This feature is evaluated for the original time series 
and for numerous realizations of a process which only ex- 
hibits the linear statistical properties of the given data. A 
significant deviation of the feature evaluated for the orig- 
inal time series from the simulated distribution suggests 
a rejection of the null hypothesis. The feature is usually 
chosen according to a specific type of alternative hypoth- 
esis on the underlying dynamics. In their Letter Palus 
and Novotna 0] chose the amplitude-frequency correla- 
tion as a property of nonlinear (driven) oscillators. But 
the rejection of the null hypothesis based on a certain 
feature does not, in general, give evidence that the spe- 
cific type of alternative that has motivated the choice of 
the feature is present. To provide evidence for a specific 
alternative one has to show that the chosen feature has 
high power to detect the violation by which it was mo- 
tivated but no power to detect other types of violations. 
Unfortunately, the null hypothesis under consideration is 
such restrictive that the possible alternatives span a huge 
class of processes, see e.g. [||,||. 

With respect to the amplitude-frequency correlation 
considered in by Palus and Novotna [Q| their Letter, for 
example, if the frequency of a second order linear stochas- 
tic process is modulated with time, the resulting pro- 
cess analytically shows an amplitude-frequency correla- 
tion g . A physically more plausible alternative hypoth- 
esis for the sun spot data arises from solar physics, see 
for review: The sun spots are an effect of the dy- 
namics of the magnetic field of the sun which exhibits a 
22 years cycle. This dynamics, a magnetohydrodynamic 
dynamo, is described by a nonlinear partial differential 
equation which is eventually stochastically driven. The 
sun spot number represent a very specific mapping from 
the spatio-temporal magnetic field to a scalar time se- 
ries. Since nonlinear driven partial differential equations 
include nonlinear driven oscillators as special cases, the 
latter can not be distinguished from the former based on 
surrogate data for the sun spots. 

Summarizing, a significant amplitude-frequency corre- 
lation is a feature of driven nonlinear oscillators, but it 
is not a specific feature of these type of processes. Thus, 
the specific alternative of a driven nonlinear oscillator can 
not be concluded from a rejection of the null hypothesis. 

Generally speaking, assuming that (1) no process in 
nature is indeed a linear, Gaussian, stationary, stochastic 


dynamical one and (2) that one is using a feature that 
is capable to detect the actual deviation, without any 
further information about the process, the only thing one 
can infer from surrogate data testing is whether there are 
enough data for the power of the test to be large enough 
to reject the by assumption untrue null hypothesis. 
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